Storage of multiple coherent microwave excitations in an electron spin ensemble 
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Electron and nuclear spins have good coherence times and an ensemble of spins is a promising 
candidate for a quantum memory. By employing holographic techniques via field gradients a single 
ensemble may be used to store many bits of information. Here we present a coherent memory 
using a pulsed magnetic field gradient, and demonstrate the storage and retrieval of up to 100 weak 
10 GHz coherent excitations in collective states of an electron spin ensemble. We further show that 
such collective excitations in the electron spin can then be stored in nuclear spin states, which offer 
coherence times in excess of seconds. 



Instead of storing information in specific locations as in 
photography and in conventional computer memory, in- 
formation can be stored in distributed collective modes, 
as in holography. Advantages include obviating the need 
for local manipulations and measurements, enhanced 
coupling to electromagnetic fields, and robustness against 
decoherence of individual members of the ensemble. This 
principle has been applied to different light-matter inter- 
faces such as atoms [IH3], ion-doped crystals [SHU]) polar 
molecules [T0UT3] , or spins [T4J [15] . Controlled reversible 
inhomogeneous broadening (CRIB) [8], or gradient echo 
memory (GEM) [7] schemes which apply external field 
gradients to address different storage modes have been 
proposed and observed in gaseous atomic samples [1] [5] 
and in ion-doped solids [6j [7] ■ 

In this Letter, we demonstrate the storage of multiple 
microwave excitations in an electron spin ensemble. The 
spin ensembles used as the storage medium are the elec- 
tron spin of nitrogen atoms in fullerene cages ( 14 N@Cgo) 
and phosphorous donors in silicon. The microwave exci- 
tations are phase encoded using a static or pulsed field 
gradient, with the latter allowing for recall in arbitrary 
order. We have stored up to 100 weak excitations in a 
spin ensemble and recalled them sequentially. We also 
demonstrate the coherent transfer of the stored multiple 
excitations between electron spin and nuclear spin, which 
will allow much longer storage times |16j . The multi- 
mode storage achieved in this way offers prospects of con- 
structing a long-lived quantum memory which could be 
used for a hybrid quantum computing architecture with 
superconducting qubits. 

A quasistatic magnetic field along the z-axis causes the 
members of the spin ensemble to process at an angular 
frequency B(z,t)fig e /H, where /x is the Bohr magneton 
and g e is the electron gyromagnetic ratio. Applying a 
magnetic gradient of strength G = dB(z, t) jdz for a time 
r consequently leads to a difference in precession angle of 
66 = {\ig e jK)Gr ■ Sz between two spins with separation 



8z along the z-axis. A gradient pulse thus maps a spin 
state with a coherent transverse magnetization (such as 
that generated by a global resonant microwave tipping 
pulse) to a spin-wave excitation associated with a wave 
number k — (fig e /h)Gr ■ z. Each further application of 
G for duration r generates a change in the wavevector of 
the global spin wave mode, by an amount k. The notion 
of one or two dimensional /c-space introduced by mag- 
netic field gradients has been widely used for a variety of 
nuclear magnetic resonance (NMR) applications [17l [18] 
including imaging |19j . 

Using this approach it is possible to store many coher- 
ences in the electron spin ensemble. If a global coherence 
is generated by a small electron spin resonance (ESR) tip- 
ping pulse 50 <C 7T, the amplitude of the coherence (given 
by the magnetization in the x-y plane) is proportional 
to sin(<5#) ~ 50, while the change in the magnetization 
along z is given by 1 — cos(50) ~ 89 2 which can be ne- 
glected. A gradient pulse converts this fcj = coherence 
into a state kf 7^ 0. These two modes are orthogonal if 
~k J-d/2 n ( z ) e ~ l< - ki ~ kj ^' Z( iz = 0, where n(z) is the density 
of spins at position z, N is the total number of spins, d 
is the extent of the sample along the z direction. This 
condition is satisfied if the total magnetization in the x-y 
plane is zero for the fc/-mode; this ensures that the kf- 
mode does not radiate as the ensemble precesses in the 
static magnetic field. Since, to first order, the total z-axis 
magnetization is unchanged by the small tipping pulse, 
a new k = coherence may be generated by a further 
small tipping pulse, and stored in another k ^ mode 
by a further gradient pulse. 

Information can be encoded on these modes by choos- 
ing the phase of the small tipping pulses, i.e. the axis 
of the tip in the x-y plane. To retrieve the information 
stored in mode fc, one needs a change in wave number 
by the amount — k to bring it back to k = mode. This 
is accomplished by applying the opposite field gradient 
— G for the same duration r. The net magnetisation of 
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this mode precesses in the external field causing it to 
radiate a signal similar to a spin echo [20] . In our ex- 
perimental setup, we generate the pulsed magnetic field 
gradient along the z-axis by placing an anti-Hclmholtz 
coil outside the resonator. The current pulses used to 
produce the field gradient pulse are highly reproducible. 
In addition to the controlled inhomogeneity of the gra- 
dient pulses, there is a random inhomogeneity caused by 
variations in the environment of the individual spins of 
the ensemble and some residual gradients in the static 
applied magnetic field. The effect of this uncontrolled 
but static inhomogeneity can be refocussed using a Hahn 
echo sequence [20]. A 7r-pulse echo also has the effect of 
inverting the wave vectors of stored collective excitations 
fc — > —k. Thus, following an odd number of refocusing tt- 
pulses, a field gradient of the original polarity G restores 
a fc-mode to the k — mode. 

Our first experimental objective is to demonstrate that 
by applying an appropriate magnetic field gradient we 
can create an orthogonal global phase mode across the 
sample. By varying the duration r (or alternatively the 
amplitude G) of the field gradient pulse, we create modes 
of varying fc, resulting in different couplings to the radia- 
tive k = mode. When we insert such a field gradient 
pulse into a Hahn echo sequence we observe a variation in 
the spin echo intensity of the N@Ceo sample, as shown in 
Fig. [T|\. A disappearance of the spin echo indicates that 
the uniform phase coherence of the spins is completely 
eliminated, and that the mode k generated by the field 
gradient pulse is orthogonal to mode k = 0. These results 
are used to obtain an appropriate duration and magni- 
tude of the gradient pulse (e.g. r = 1.3 fis and ~ 30 
mT/m) that will be used in the subsequent experiments 
for storing multiple microwave excitations. 

The elimination of phase coherence by using magnetic 
field gradients has been exploited in NMR experiments 
[2"Tl |2"2"] where the field gradient was used as a phase co- 
herence eraser. Here, however, we use the pulsed field 
gradient to encode and label excitations stored in the 
multi-mode memory. In Fig. [Tj3 we show that a sin- 
gle microwave excitation generated by a microwave tt/2 
pulse and stored into a fc-mode orthogonal to the fc = 
mode, can be retrieved by applying the same field gra- 
dient pulse after the refocusing 7r-pulse. The recovered 
spin echo appears in the same place as for a standard 
Hahn echo sequence, with a fidelity of 92%. 

A key resource for both classical and quantum com- 
putation is a memory allowing reading and writing of 
registers. The memory is particularly useful if regis- 
ters are accessible in arbitrary order. In our scheme, we 
achieve this by converting the register to be accessed to 
the fc = mode using appropriate gradient pulses. The 
gradient pulses also shift all other registers in fc-space by 
the same amount; the memory can be thought of as a 
"Turing tape" residing in fc-space. 

Figure [2] shows an experiment in which we store two 




FIG. 1: Single microwave pulse stored in a mode or- 
thogonal to fc = mode across a spin ensemble using 
a field gradient pulse (A): The integrated intensity of the 
spin echo is plotted against the integral of the amplitude of 
the gradient pulse over time. Considering the profile of the 
sample, the fitted curve is of the form Bi(kro)/kro (see Meth- 
ods), corrected with a second order radial inhomogeneity of 
the microwave field. Points of zero echo intensity correspond 
to modes orthogonal to k — mode. (B): A second identical 
gradient pulse is applied after the refocussing it pulse to shift 
the stored pulse back iinto k — mode, which then radiates 
to give an echo. Sample: 14 N@C6o at room temperature. 

weak microwave excitations of different phase (Pi and 
P2) in the spin ensemble and recall them in either order. 
Each excitation is followed by an encoding field gradient 
pulse G. Since the field gradient pulses act on all the 
microwave excitations stored in the ensemble, the first 
excitation Pi is effectively stored into the mode 2fc(G, r) 
while the second P2 is stored in mode k(G, r). Therefore, 
after the first refocusing 7r-pulse, we can restore and refo- 
cus either Pi or P 2 to the fc = mode by applying either 
two or one gradient pulses respectively. This encoding 
method can be readily generalized to the case of multi- 
ple microwave excitations: by using pulsed field gradients 
we can address any of the stored excitations and retrieve 
them individually. 

In the limit of single quantum excitations stored in or- 
thogonal fc-modes, there is no crosstalk between stored 
registers. However, in our experiment the excitations 
that we store and retrieve are far from single quantum ex- 
citations. Therefore, we investigate the crosstalk between 
two stored microwave excitations as a function of their 
intensities. We compare the echo intensity I of one recov- 
ered excitation in the presence of another excitation with 
the echo intensity Iq in the absence of another excitation, 
and calculated the fractional change D = (Iq — I) /Jo- 
Supposing that the two microwave excitations are stored 
in two orthogonal modes, we would expect a decrease in 
the echo of the second excitation (Pj) when increasing 
the intensity of the first excitation (Pi) because the net 
magnetization along the z direction is partly consumed 
by Pi, leading to D2 = 1 — cos(#i) where 6\ is the tip- 
ping angle for Pi. We would also expect a decrease in the 
echo of Pi with the increase of the intensity of P2 due to 
the partial refocusing effect of P2, which is a result of the 
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FIG. 2: Recalling two stored microwave pulses in an 
arbitrary order using pulsed field gradient (A): Two mi- 
crowave pulses recalled in the same order as they are stored. 
(B): Two pulses recalled in the inverse order. Transients are 
shown respectively for Pi and P2 applied in +x and —x, +x 
and +y direction. In each of the four diagrams we show the 
real part (+x direction, black) and imaginary part (+y di- 
rection, red) of the signal. As the refocusing pulse is in y 
direction, the imaginary part of the signal between the first 
and second refocusing pulse (orange) is inverted to keep the 
echo and free induction decay (FID) having the same sign. 
In this experiment Pi and P2 are n/6 pulses. They are kept 
relatively strong so that we can observe a clear echo. Sample: 
14 N@C6o at room temperature. 



fact that any two tipping pulses generate some refocusing 
(even if the second pulse is not a 7r-pulse as required for 
optimal refocusing). This leads to D\ = [1 — cos(#2)]/2, 
where #2 is the tipping angle for P 2 . We found that Di 
and Z?2 follow the expected dependencies within 0.5%, for 
0\ and 62 each varying independently in the range to 
7r/2, confirming that we have identified the key crosstalk 
mechanisms (see Supplementary Information). Since the 
error D depends on the tip angles to second order for 
both mechanisms, it quickly becomes negligible for small 
tip angles. 

Apart from this crosstalk, errors in the echoes can also 
arise from errors in the microwave pulses [23 . An inac- 
curate refocusing 7r pulse will introduce extra excitations 
in the spin ensemble and pollute the stored information, 
but there are well established methods for mitigating the 
dominant instrumental imperfections |24j . Alternatively, 
such errors could be avoided if we were able to invert the 
inhomogeneity without using microwave pulses. Tech- 
niques of this kind have been developed for multi-mode 
memories with atomic ensembles, e.g. CRIB and atomic 
frequency comb |25j . 

Having demonstrated that we are able to store multiple 
excitations and that the crosstalk can be kept small for 
weak excitations, we address the question of how many 
excitations it is practically possible to store. In this ex- 
periment we apply a static field gradient to an ensemble 
of P-donor spins in silicon by placing a small permanent 
magnet close to the resonator, and apply a train of weak 
(less than 0.0l7r) excitations at intervals in time such that 
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FIG. 3: Multiple microwave pulses stored in multi- 
modes of the ensemble using constant field gradient 

(A): Storage of 100 weak microwave pulses. Each pulse has a 
flip angle of less than 0.01 tt and is applied arbitrarily in the 
+x or —x direction. The FID of the sample is shorter than 2 
/is in this constant field gradient, so the microwave pulses are 
separated by 3 fis. After the refocusing pulse the stored pulses 
are recovered in reverse time order. (B): 8 weak microwave 
pulses are stored in the electron spin ensemble using constant 
gradient. All the different k modes are transferred into the 
nuclear spin, after applying an rf 7r-pulse, transferred back to 
the electron spin. The transient shows the 8 echoes (enclosed 
in the red dashed box in the pulse sequence) observed corre- 
sponding the 8 stored excitations. Sample: P-doped 28 Si at 
9 K. 



each excitation becomes orthogonal to the k = mode 
before the application of the subsequent one, i.e. the time 
separation of pulses exceeds the ensemble dephasing time 
T£ . The microwave pulses are much shorter than Tjf, 
ensuring that the bandwidth is sufficient to excite the 
whole ensemble uniformly. Fig. [3]^ shows the storage and 
retrieval of the 100 weak microwave pulses individually. 
The pulses are applied arbitrarily in +x and — x direction 
to represent a register of information. As the field gradi- 
ent is constant, only one refocusing 7r-pulse is sufficient 
to recall all the stored excitations. The wave number reg- 
ister forms a quantum memory stack, and the excitations 
are recalled in reverse order of there storage [17] . 

The memory time of the information is limited by the 
T 2 decoherence time of the electron spin (about 450 ^s 
for P-donor spin at 9K), as seen in the decrease of the 
intensity of the echoes in Fig. [3]A_. On the other hand, 
the shortest separation between each microwave excita- 
tion is the T£ in the presence of the gradient (in this case 
2 /zs). These two figures set the limit for the number 
of microwave excitations that can be stored in the en- 
semble. Although a memory time of several hundreds of 
microseconds is already much longer than the best super- 
conducting qubits, and times as long tens to hundreds of 
milliseconds have been reported [26-428] . the stored ex- 
citations can be transferred to nuclear spin in order to 
benefit from even longer nuclear spin coherence times in 
excess of seconds \Tjg\ . Fig. [3(3 shows the storage of eight 
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weak microwave pulses using a constant field gradient. 
We then transfer the entire state of the electron spin 
ensemble (i.e. all the microwave excitations stored in dif- 
ferent k modes) simultaneously to a nuclear spin degree 
of freedom. That this transfer is a transition between hy- 
pcrfme states of each individual atom and it can hence be 
carried out on a time scale given by the hypcrfinc inter- 
action strength, which is in the order of microseconds in 
our system. As the excitations are now in nuclear spins, 
we achieve the refocusing with a radio frequency 7r-pulse 
resonant with the nuclei instead of a microwave 7r-pulse. 
We then transfer the coherences back to the electron spin 
and echoes of the eight original microwave excitations are 
observed. 

Our experiments with small tipping angles demon- 
strate how a single electron spin ensemble could be 
used as a multi-mode solid-state memory. The pre- 
requisite techniques for extending this work to the 
quantum regime, synthesizing arbitrary single photon 
states [53] and non-destructively detecting single cav- 
ity photons [3D], have now been achieved in experiments 
with superconducting qubits. Recent experiments have 
also demonstrated that large couplings can be achieved 
between spin ensembles and cavities [311 132] ■ Further- 
studies are necessary to show that the demonstrated co- 
herence of the qubits and spins can be maintained in a 
combined system. 

In summary, we have demonstrated holographic stor- 
age of microwave pulses in an electron spin ensemble. 
By using magnetic field gradients the microwave pulses 
are encoded in different collective modes of the spin en- 
semble, and are retrieved individually. More robust stor- 
age can be achieved by transferring the coherence from 
electron spin to a coupled nuclear spin. These results 
show the prospect of using spin ensembles as a memory 
medium and implementing a quantum computing scheme 
with hybrid physical systems. 
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Supplementary Information 

In Fig. 4 we show the crosstalk between different modes 
A: as a function of the intensity of the stored microwave 
excitations. We plot both experimental data and theo- 
retical calculations. The curves for the theoretical results 
are D\ = (1 — cos #2)/2, -D2 = 1 — cos 6\, with Q\ and O2 
the flip angles of Pi and P2. The deviation between the 
experimental data and the theoretical curves is less than 
0.5%. 




FIG. 4: Difference in the recovered echo intensity (D) 
of one microwave excitation caused by another excita- 
tion (A): D for the echo of the first microwave pulse (Pi) Ei 
as a function of the flip angle of the second microwave pulse 
P2. (B): D for E2 as a function of the flip angle of Pi. In 
both cases the durations of Pi and P2 are swept from less than 
0.027T to 0.6-7T, however D was independent of the intensity of 
the microwave pulse whose echo intensity was measured. 



